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SUMMARY 


Modern  numerical  methods  for  the  optimization  of  large  dis¬ 
cretized  systems  are  now  well  developed  and  highly  efficient  in 
the  case  of  thin  vailed  elastic  structures  m  c  a  e  1  e  d  by  finite* 
elements.  However,  this  is  not  yet  true  for  structures  whose 
components  are  subject  simultaneously  to  bending  and  extension 
loads.  In  this  paper,  the  idea  of  Generalized  Optimality  Cri¬ 
terion  (GOC),  set  forth  in  previous  final  scientific  reports 
for  bar,  membrane  and  pure  bending  elements,  is  extended  to  deal 
with  general  beam  and  flat  shell  elements.  The  modifications 
brought  to  the  GOC  result  in  explicit  approximations  for  the 
behavior  constraints  that  are  still  correct  up  to  the  first,  orde 
but  that  exhibit  a  more  complex  algebraic  form.  Indeed  these 
explicit  expressions  are  no  longer  merely  linear  in  the  recipro, 
design  variables.  However  they  continue  to  be  addi lively  sepa¬ 
rable  and  therefore,  dual  methods  remain  fully  applicable,  just 
as  in  the  original  statement  of  the  GOC  approach.  Numerical 
examples  will  be  offered  to  demonstrate  the  efficiency  of  the 
method  presented. 


INTRODUCTION 


1  . 


The  optimum  design  of  any  significant  structure  is  the  result 
of  a  delicate  compromise  between  many  complex  factors.  Some  are 
rational  and  can  be  quantified,  such  as  the  strength  of  the 
structure,  its  natural  frequencies,  its  weight,  ...  Some  are 
just  as  rational  but  are  difficult  to  quantify,  such  as  the  ex* 
perience  of  the  designer  in  a  given  technology.  Some  others  arc 
much  less  rational,  like  styling,  but  are  just  as  important  for 
the  final  goal  of  the  process,  which  is  the  marketing.  Naturally 
a  good  designer  considers  structural  optimization  as  a  technique 
that  should  take  into  account  all  possible  aspects  of  the  design. 
Consequently  the  designers  are  often  reluctant  to  the  concepts 
of  structural  optimization  developed  in  connection  with  finite 
element  programs. 

However  a  more  detailed  examination  of  the  design  process 
leads  to  isolate  a  phase  that  appears  frequently,  during  which 
the  shape  of  the  structure  is  more  or  less  frozen  and  the  problem 
is  limited  to  giving  adequate  dimensions  to  the  various  members. 
Such  a  situation  is  often  encountered  in  the  aerospace,  naval 
or  automobile  industries,  where  the  external  shapes  are,  to  a 
large  extent,  dictated  by  aero-  or  hydrodynamic  considerations, 
or  by  styling,  while  internal  forms  are  often  determinated  by 
various  other  non  structural  requirements.  If  the  ultimate  goal 
of  the  designer  can  be  identified  as  corresponding  to  the  minimi¬ 
zation  of  an  explicit  function  of  the  member  sizes,  and  if  the 
limitations  on  the  design  can  be  defined  as,  eventually  implicit, 
functions  of  the  member  sizes  too,  such  as  displacements,  stresses, 
eigenfrequencies,  etc...,  then  the  problem  is  tractable  by  auto¬ 
matic  algorithms.  They  allow  the  designer  to  speed  up  signifi¬ 
cantly  this  part  of  the  design  process  and  to  explore  more  sys¬ 
tematically  the  various  feasible  designs. 

The  optimization  problems,  which  in  fact  should  be  called 
automatic  sizing  problems,  are  especially  crucial  when  complex 
structural  forms  are  involved,  when  flexural  forces  cannot  be 
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neglected,  and  when  composite  materials  such  as  reinforced 
resins  are  employed.  In  these  cases  it  becomes  difficult,  if 
not  impossible,  for  the  designer  to  have  an  intuitive  under¬ 
standing  of  the  structural  mechanics  that  is  sufficient  to  lead 
to  op t  i  mal  s i z i ng  of  the  various  members.  Furthermore,  the  de¬ 
signer  is  most  of  the  time  unable  to  take  into  account  global 
constraints  in  the  structure,  like  global  flexibility,  restric¬ 
tion  on  displacements,  frequencies  of  vibration,  global  buckling 
modes,  etc...  It  is  only  possible  to  verify  a  posteriori  that 
such  constraints  are  satisfied.  Again  these  global  constraints 
become  more  important  in  the  context  of  highly,  indeterminate 
structures.  In  the  aerospace  industry,  the  necessity  of  desi¬ 
gning  high  performance  structures  has  motivated  significant  re¬ 
search  efforts  to  derive  algorithms  permitting  a  rapid  and  sys¬ 
tematic  exploration  of  the  design  space  to  determine  the  optimum 
material  utilization. 

It  is  worth  pointing  out  that  optimization  methods  should 
be  considered  as  especially  useful  in  the  preliminary  design 
phase.  Using  them  when  the  design  is  practically  frozen,  with 
the  hope  of  an  ultimate  improvement,  is  often  disappointing. 

This  is  due  to  the  fact  that  the  optimization  of  a  detailed  design 
implies  the  formulation  of  a  large  number  of  constraints,  some 
of  which  are  not  easily  quantified.  At  the  preliminary  design 
stage,  however,  the  constraints  are  usually  more  global  and 
therefore  more  easily  handled  by  the  available  formulations. 

The  structural  optimization  problem  considered  in  this  report 
consists  of  the  weight  minimization  of  a  finite  element  model 
with  fixed  geometry  and  material  properties.  The  design  variables 
are  taken  as  the  transverse  sizes  of  the  structural  members, 
namely,  the  cross-sectional  areas  of  bar  and  beam  elements  and 
the  thicknesses  of  membrane,  plate  and  flat  shell  elements. 

The  mathematical  programming  problem  to  be  solved  has  the  fol¬ 
lowing  form  : 
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The  a^?s  denote  the  n  design  variables.  They  correspond  to 
member  sizes  of  either  individual  finite  elements,  or,  if 
design  variable  Linking  is  used,  of  groups  of  finite  elements. 
The  structural  weight  W  is  a  linear  objective  function,  because 
the  are  constant  coefficients  representing  the  specific 

weight  times  either  the  element  length  (bars  and  beams)  or  the 
element  area  (membranes,  shear  panels  ,  plates  and  flat  shells). 
The  inequalities  (2)  are  the  behavior  constraints,  which  impose 
limitations  on  quantities  describing  the  structural  response, 
for  example,  the  stresses  and  displacements  under  multiple 
loading  cases,  the  natural  frequencies,  the  buckling  loads, 
etc...  The  design  variables  are  also  subjected  to  the  side 
constraints  (3),  where  and  are  lower  and  upper  limits  that 
reflect  fabrication  and  analysis  validity  considerations. 

Standard  minimization  techniques  can  be  applied  to  the  non¬ 
linear  programming  problem  (1-3).  However  this  problem  exhibits 
some  characteristics  that  make  it  complicated  when  practical 
structural  design  applications  are  considered.  The  essential 
difficulty  arises  from  the  implicit  nature  of  the  behavior 
constraints  (2),  in  that  their  precise  numerical  evaluation  for 
each  particular  design  requires  a  complete  finite  element  ana¬ 
lysis.  Since  the  solution  scheme  is  iterative,  it  involves  a 
large  number  of  structural  reanalyses.  Therefore  the  computa¬ 
tional  cost  often  becomes  prohibitive  when  large  structural 
systems  are  dealt  with.  However  a  powerful  design  procedure 
has  now  emerged,  which  consists  in  replacing  the  initial  problem 
with  a  sequence  of  simple  explicit  problems.  In  the  next  sec¬ 
tion  this  approach  will  be  briefly  reviewed  by  restricting  the 
formulation  to  thin  walled  structures  idealized  by  bar/membrane 
elements.  A  much  more  detailed  presentation  can  bo  found  in  a 
previous  report  [  1  ]  .  It  will  be  shown  that  the  behavior  cons¬ 
traints  can  be  approximated  either  by  using  virtual  load  consi¬ 
derations  (optimality  criteria  approach)  or  by  using  first  order 
Taylor  series  expansion  with  respect  to  the  reciprocal  design 
variables  (mathematical  programming  approach).  Applying  a  dual 
solution  scheme  to  each  explicit  problem  generated  in  sequence 
naturally  introduces  the  concept  of  a  generalized  optimality 


criterion. 


4  . 


Subsequently  sections  3  and  4  wiLL  be  concerned  with  struc¬ 
tural  models  that  are  capable  of  carrying  flexural  forces.  For 
beams  and  plates  in  pure  bending,  adequate  intermediate  varia¬ 
bles  can  be  selected,  in  terms  of  which  high  quality  explicit 
approximations  for  the  behavior  constraints  can  still  be  gene¬ 
rated  by  linearization.  The  idea  of  generalized  optimality 
criterion  remains  fully  valid  and  it  keeps  its  interpretation 
in  terms  of  energy  densities  in  the  structural  members.  Sect i on 
3  is  a  summary  of  results  presented  in  a  previous  report  (ref. 

[2  1,  section  6). 

The  essential  problem  is  addressed  in  section  4.  It  consists 
of  the  establishment  of  the  generalized  optimality  criterion 
approach  in  the  general  case  where  the  structural  members  work 
both  in  extension  and  flexion  (beam  and  flat  shell  elements). 

For  displacement  constraints,  it  is  no  longer  possible  to  select 
a  suitable  intermediate  variable  for  the  linearization  process, 
however,  the  virtual  load  procedure  permits  obtention  of  high 
quality,  first  order  explicit  approximations  of  the  behavior 
constraints.  It  will  also  be  shown  how  to  proceed  for  stress, 
frequency  and  buckling  constraints. 

The  explicit  approximations  still  exhibit  a  separable  alge¬ 
braic  form  and  therefore  dual  methods  remain  applicable.  However 
some  difficulties  might  happen  due  to  the  lack  of  convexity  in 
the  explicit  sub  problem.  Section  3  will  describe  how  to  solve 
the  explicit  sub  problems  by  using  dual  methods.  Finally  seine 
applications  will  be  offered  in  section  6  to  illustrate  the  power 
and  generality  of  the  approach  presented. 


2.  GENERALIZED  OPTIMALITY  CRITERION  FOR  THIN  WALLED  STRUCTURES 


This  section  summarizes  some  results  obtained  in  a  previous 
work  (  1  ]  for  structural  models  made  up  of  bar,  membrane  and 

shear  panel  elements,  which  are  quite  adequate  for  idealizing 
thin  walled  structures  subjected  mainly  to  extension  loading. 
For  this  class  of  finite  element  models,  the  structural  stif¬ 
fness  matrix  exhibits  a  linear  form  in  the  design  variables  : 
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.  Z  , 

x  =  i 


a.  K  . 
1  1 


(  •'  ) 


w he  re  K^,  a  matrix  of  constants,  represents  the  stiffness  matri 
of  the  ith  element  when  =.  1.  For  simplicity,  the  following, 

discussion  is  restricted  to  problems  involving  constraints  on 
static  stresses  and  displacements,  in  which  case  the  behavior 
constraints  (2)  can  be  written 


h.(a)  ;  u  .  -  u.(a)  >  0 

J  J  j 


(  r)  ) 


where  u.  denotes  an  upper  bound  to  a  response  quantity  Uj(n) 
(stress,  nodal  displacement,  relative  displacement). 


Most  of  the  optimality  criteria  approaches  (  e  .  g  .  |  3  ]  ),  as 
well  as  the  generalized  optimality  criterion  (GOC)  set  forth  in 
Ret.  [  1  ] ,  use  the  virtual  load  technique  to  generate  explicit 

approximations  of  the  stress  and  displacement  constraints.  In¬ 
troducing  a  virtual  load  vector  conjugated  to  the  response  quan¬ 
tity  Uj  (unit  load  for  a  nodal  displacement),  it  follows  that 
u j  can  be  expressed  as  the  sum  of  the  contributions  of  each  fi¬ 
nite  element  : 
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Ln  these  expressions  q  and  are  respectively  the  re,!  1  and 

virtual  displacement  vectors  and  K.  is  the  clement  stiffness 

matrix  appearing  in  ( 4 )  .  It  can  be  seen  from  (7)  that  t  h  e  <  ■  <  >  e  t - 

f  i  c  i  o  ills  c  .  .  a  r  e  r  e  ]  a  ted  to  the  virtual  s  t  r  a  i  n  e  u  c  r  y  v  densities 
ij 

in  the  structural  members.  The  c  .  .  1  s  are  constant  coefficients 

i  J 

in  the  case  of  a  statically  determinate  structure,  so  that  (f>) 
represents  then  the  exact  explicit  form  of  the  response  quan- 
citv  u j .  In  the  case  of  a  statically  indeterminate  structure, 
the  c^jrs  depend  implicitly  on  the  design  variables,  because 
structural  redundancy  produces  redistribution  of  the*  internal 
forces  when  the  member  sizes  are  modified.  Therefore  the  inf¬ 
lowing  explicit  constraints 
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constitute  in  general  approximate  forms  of  the  original  cons¬ 
traints  (5).  As  shown  in  Ref.  |  1  ],  the  basic  idea  in  the  op¬ 
timality  criteria  approach  can  be  viewed  as  transforming  the 
initial  implicit  problem  into  a  sequence  of  explicit  subproblem 
Each  explicit  problem  results  from  replacing  the  behavior  cons¬ 
traints  (2)  by  their  approximate  forms  (8). 


On  the  other  hand,  the  mathematical  programming  approach  to 
structural  optimization,  after  a  period  of  unefficiency,  has 
finally  evolved  into  a  powerful  and  now  well  established  design 
procedure  which  is  also  based  upon  explicit  approximations  ot 
the  behavior  constraints  [4,  5,  6  J  .  The  key  idea  is  to  linea¬ 

rize  the  behavior  constraints  with  respect  to  the  reciprocal 
design  variables 


Justification  for  this  change  of  variables  lies  in  the  fact 
that  the  constraint  surfaces  can  be  shown  to  be  very  shallow 
and  close  to  planes  in  the  reciprocal  design  variable  space. 
Therefore  the  linearized  forms  of  the  constraints  are  usually 
high  quality  approximations.  They  are  obtained  by  using  first 
order  Taylor  series  expansion  in  terms  of  the  reciprocal  va¬ 
riables  x.  : 
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where  the  uppcrscript  °  denotes  quantities  evaluated  at  the 
actual  design  point  x°,  w  tie  re  the  structural  analysis  is  p  •  ■  re¬ 
formed.  Nc;e  that  the  finite  element  analysis  capability  m  n  n  t 
include  auxiliary  sensitivity  analyses  for  evaluating  the  lirsi 
partial  derivatives  of  the  response  quantities.  Most  often  the 
well  known  pseudo-loads  technique  is  employed  |  7  ]  . 


It  has  been  shown  in  a  previous  report  |  1  | 

approximations  of  the  behavior  constraints  used 
optimality  criteria  and  mathematical  programming 
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1 0 ,  respectively)  are  identical.  Indeed  the 
densities  c  ±  j  employed  in  the  optimality 
are  nothing  else  that  the  gradients  of  the 
with  respect  to  the  reciprocal  variables 
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Furthe rmo  re  the  de  f i n i t i on  of  the  c .  .  ' s  following  from  the  v  i  r  - 

i  J 

t  u  a  1  load  technique  (see  E  q .  7)  clearly  indicates  that 
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Therefore  (10)  can  be  rewritten 

fr.(x)u.  -.Z.c?.  x.  >0  (l<) 

j  j  i  =  i  i  j  i 

which  is  equivalent  to  (8)  when  recast  in  terms  of  the  direct 
variables  a^.  It  is  thus  apparent  that  a  unified  approach  to 
structural  weight  minimization  of  finite  element  systems  has 
emerged,  which  consists  in  replacing  the  initial  problem  (1-3) 
with  a  sequence  of  explicit  approximate  -  or  linearized  -  pro¬ 
blems  of  the  following  form  : 
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The  GOC  statement  results  from  writing  the  KUHN -TUCKER  optimal i tv 
conditions  for  the  problem  (14-16).  This  yields  an  explicit 

expression  for  the  design  variables  a.  in  terms  o  I  the  1  a  g  r  a  n  ;■  i  .1  :i 

multipliers  r.  associated  with  the  behavior  curis  t  mints  (  I '»  ) 

(see  Ref  |  6  ,  H  |  for  more  details) 

I  m  1  / 1  y  m  _  ) 

a  .  =  ( -- —  .  T.  c  .  .  r  .  )  ~  if  Z  .  a  T  <  .2  c  .  .  r  .  •  ■  .  a  T  (  1  7  ) 
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Of  course  the  lagrangian  multipliers  must  be  uonnegative,  more 
precisely  >  they  must  satisfy  the  complementary  conditions  : 
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In  order  to  compute  the  lagrangian  multipliers  satisfying  (20— 
21),  an  interesting  approach  is  to  resort  to  dual  methods,  which 
leads  to  maximize  the  lagrangian  function  considered  as  a  func¬ 
tion  of  the  lagrangian  multipliers  -  or  dual  variables  -  only 
(see  Refs  (  1,  6,  8  ]  and  section  5).  Once  the  solution  of  the 

dual  problem  has  been  found,  the  corresponding  optima]  design 
variables  are  easily  computed  from  the  explicit  optimality  cri¬ 
teria  relations  (17-19).  Note  that  the  design  variables  can  be 
separated  into  a  group  of  active  (or  free)  variables  (see  E  q .  17) 
and  a  group  of  passive  (or  fixed)  variables  (see  E  q  s  .  18  and  19). 
This  subdivision  into  active  and  passive  design  variable  groups 
is  classical  in  the  optimality  criteria  approaches  |  3  ,  9  ].  It 

corresponds  to  the  fact  that  the  dual  space  -  i.e.  the  space  of 
the  lagrangian  multipliers  r^  -  is  partitioned  into  several  re¬ 
gions  separated  by  planes  across  which  the  second  derivatives 
of  the  dual  function  are  discontinuous  [6,  8  1  . 


The  whole  process  of  combining  the  linearization  of  the 
behavior  constraints  with  respect  to  the  reciprocal  design 
variables  and  a  dual  solution  scheme  can  be  viewed  as  a  gene¬ 
ralization  of  the  optimality  criteria  approaches.  It  is  impor¬ 
tant  to  mention  that  this  basic  approach  of  converting  the 
initial  problem  into  a  sequence  of  explicit  sub  problems  is  now 
widely  recognized  [6,  10  ]  and  it  is  routinely  employed  for 

large  scale  industrial  applications  [  11  ]  . 


'3.  purl:  bunding  elements 


In  this  section  attention  is  focused  on  discretized  models 
made  up  of  pure  beam  and  plate  elements  subjected  to  flexural 
loads  only  (for  more  details,  see  section  6  of  Ref.  |2  ]  )  . 

The  stiffness  matrix  of  such  a  bending  element  is  usually  not 
merely  proportional  to  its  cross-sectional  size  and  therefore 
the  optimization  strategy  reviewed  in  the  previous  section  must 
be  modified.  The  way  to  deal  with  a  beam  element  subjected  to 
uniaxial  bending  depends  upon  the  relationship  between  the  prin¬ 
cipal  moment  of  inertia  I  and  the  cross-sectional  area  a.  A 
wide  variety  of  situations  is  taken  into  consideration  by  adop¬ 
ting  the  following  relation  : 

I  =  C  aP  (22) 

where  c  is  a  constant  that  depends  only  on  the  shape  of  the 
beam  cross-section  and  p  is  a  positive  number. 

Most  of  the  time  p  is  taken  as  an  integer  number,  equal  in 
1,  2  or  3.  The  case  p  =  1  corresponds  to  thin  walled  beams, 
for  examp  le ,  sandwich  beams,  pipes  with  fixed  diameter  and  va¬ 
riable  thickness,  etc...  The  GOC  approach  of  section  2  remains 
then  fully  applicable,  since  the  stiffness  matrix  continues  to 
be  linear  in  the  design  variables.  The  case  p  =  2  is  that  of 
beams  with  uniformly  varying  cross-section.  The  shape  of  the 
cross-section  is  kept  unvariant  while  its  area  is  modified  during 
redesign  (dilatation  or  contraction).  Finally  the  case  p  =  3 
is  concerned  with  beams  having  full  cross-section  whose  height 
varies  while  other  sizes  are  fixed.  For  a  beam  subjected  to 
pure  bending,  the  flexural  rigidity  is  proportional  to  the  mo¬ 
ment  of  inertia  and  therefore,  in  a  finite  element  context,  the 
structural  stiffness  matrix  exhibits  the  foil  ow  i  up,  explicit 
form  in  terms  of  the  cross-sectional  areas  : 

k  -  *,  -  J,  K.  ,  ,,  >  o  <r«) 

where  each  matrix  K  ^  is  independent  of  the  design  variables 


a  .  . 
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With  regard  to  plate  elements  subjected  to  pure  bending, 
two  cases  must  be  disting uished.  The  first  cast'  is  that  of 
sandwich  plates  with  constant  core  thickness.  The  sheet  thick- 
nesses  constitute  then  the  design  variables.  Consequently  tin* 
stiffness  matrix  continues  to  depend  linearly  on  the  design 
variables  and  the  GOC  relations  (17-21)  remain  fully  applicable. 
The  second  case  is  concerned  with  full  plates  with  variable 
thickness.  The  stiffness  is  then  proportional  to  the  cube  ol 
the  thickness  and,  in  an  assembling  of  plate  elements,  the  re¬ 
lation  (22)  must  be  chosen  with  p  =  3. 

By  assuming  that  the  structural  discretization  is  made  up 
entirely  of  elements  of  the  same  type,  the  stiffness  matrix 
exhibits  the  form  (23),  where  p  takes  on  the  same  value  foi 
each  member.  In  these  circumstances,  the  GOC  can  be  derived 
just  as  in  the  case  of  thin  walled  structures,  by  adopting  a 
change  of  variables  tending  to  reduce  the  nonlinear  character 
of  the  constraints  : 


a  . 

l 

The  next  step  is  to  linearize  the  constraints  with  respect  to 
the  new  variables  x^,  which  requires  gradient  evaluation  (see 
bq.  10).  Restricting  again  the  discussion  to  stress  and  dis¬ 
placement  constraints,  it  is  easily  shown  (see  Ref.  |2  |)  that 
the  first  order  Taylor  series  expansion  (10)  reduces  to  the 
form  (13),  or,  when  written  in  terms  of  the  direct  variables  : 
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c^j  coefficients  can  be  interpreted  as  the  gradients  of  the 


response  quantities  with  respect  to  the  intermediate  variables 

defined  in  (24),  but  they  can  also  be  related  to  the  virtual 

strain  energies  e..  in  the  structural  members  : 
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where,  by  definition. 


f  :> . 
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In  this  connection  it  should  b e  recognized  that  virtual  load 
considerations  could  directly  be  e  mp 1  o  y  e  d  to  derive  the  explicit 
approximations  (25),  instead  of  resorting  to  first  order  Taylor 
series  (see  next  section). 


The  optimality  criteria  equations  are  very  similar  to  ( l / — 
19),  provided  that  care  is  taken  of  the  exponent  p  appearing, 
in  (25).  The  basic  redesign  relations  for  the  active  design 
variables  must  read  as  follows  : 
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(28) 


where  it  is  understood  that  the  dual  variables  r .  (i.e.  the 

lagrangian  multipliers)  must  satisfy  the  complementarity  con¬ 
ditions  (see  Eqs .  20  and  21).  A  physical  interpretation  of  the 

optimality  criterion  is  obtained  by  introducing  the  virtual 
strain  energy  densities  per  unit  weight  : 
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where  e  .  .  is  defined  in  (27).  In  terms  of  the  ?  s  ,  the  o  p  t  i 

ij  i.) 

mality  criterion  (28)  takes  the  "energetic"  form 
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In  the  special  case  where  only  one  displacement  constraint 
is  specified,  the  optimality  criterion  states  that  the  virtual 
strain  energy  density  must  be  the  same  in  each  element.  In 
this  simple  case,  it  is  possible  to  solve  analytically  the  ex¬ 
plicit  problem  and  to  derive  explicit  redesign  relations  in  terms 
of  known  quantities.  The  active  design  variables  can  be  shown 
to  be  given  by 
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while  the  remaining  passive  variables  are  fixed  to  an  upper  or 
a  lower  limit  (uq  denotes  the  contribution  of  these*  passive 
variables  to  the  displacement  constraint  u  ^  u).  It  is  worth 
mentioning  that  (31)  is  well  suited  for  the  design  of  p lores 
in  bending  with  a  single  displacement  constraint.  Since  then 
p  =  3  the  redesign  relation  (31)  involves  the  fourth  root  of  the 
coefficients  c^,  rather  than  the  third  root  as  employed  in 
Ref.  [12  ]  on  an  intuitive  basis.  Note  also  that  by  taking 
p  =  1  in  (31),  conventional  redesign  relations  are  recovered, 

which  were  devised  for  trusses  [3  ]  ,  sandwich  beams  (  1  3  1  ,  *  ■  t c .  .  . 

It  can  be  concluded  that  the  generalized  optimality  criteria 
approach  can  easily  be  extended  to  deal  with  pure  bending  ele¬ 
ments  by  defining  adequate  intermediate  variables.  High  quality 
explicit  approximation';  of  the  behavior  constraints  can  sti  1  1  be 
generated  an^  the  resulting  GOC  keeps  its  interpretation  in  terms 
of  energy  densities  ir  the  structural  members.  Seeking  the  de¬ 
sign  variables  that  satisfy  the  GOC  at  each  redesign  stage  can 
still  be  achieved  efficiently  by  resorting  to  dual  methods,  be¬ 
cause  the  explicit  approximate  problem  remains  separable  and 
strictly  convex  when  expressed  in  the  intermediate  design  vari¬ 
ables  x  .  . 

l 


4 .  FLEX tON-EXTENS ION  ELEMENTS 


When  flexion  and  extension  loadings  act  simultaneously  with 
comparable  intensity  at  the  element  level,  the  definition  (  2  3 ) 
of  the  stiffness  matrix  can  no  longer  characterize  the  structu¬ 
ral  model  with  sufficient  accuracy.  To  help  fix  ideas,  consi¬ 
der  a  flat  shell  element  made  up  of  a  membrane  and  a  plate  sta¬ 
cked  together.  The  stiffness  matrix  of  such  a  11  at  shell  ele¬ 
ment  exhibits  the  form 


K  . 
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(rz) 


where  K  f  ^  and  K  f  are  constant  matrices.  As  a  result,  in 
l  i 

the  GOC  approach,  if  the  constraints  are  linearized  with  respect 
to  the  reciprocal  design  variables  (9),  their  first  order  expli¬ 
cit  approximations,  given  by  expressions  similar  to  (8),  will 
be  of  high  quality  only  if  the  structural  members  behave  main  1 v 
in  extension.  On  the  other  hand,  if  the  bending  behavior  is 
dominant,  it  is  better  to  adopt  the  change  of  variables  (24), 
yielding  first  order  explicit  approximations  of  the  form  (2r>) 
(with  p  =  3).  As  a  matter  of  fact,  the  true  situation  is  usu¬ 
ally  a  combination  of  extension  and  bending.  In  a  practical 
structure,  some  members  work  mainly  in  extension,  some  in  flexion 
and  others,  simultaneously  in  flexion  and  extension.  Whence 
the  idea  of  using  the  following  explicit  approximations,  which 
should  be  valid  in  any  situation  : 
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where  the  coefficients  cfP  and  are  considered  as  constant 

i  J  1.1 

throughout  the  redesign  phase. 


Because  it  is  no  longer  possible  to  select  appropriate  in¬ 
termediate  design  variables,  the  explicit  approximations  (VO 
cannot  be  obtained  by  merely  using  first  order  Tavlor  series 
as  in  the  case  of  pure  bending  elements.  However  an  essential 
requirement  is  that  these  approximations  remain  correct  up  to 
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Che  first  order,  despite  the  fact  that  they  do  not  result  from 
a  strict  linearization  process.  In  other  words  the  following 
equality  must  hold  ; 
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This  condition  insures  that,  at  the  optimum,  the  solution  to 
the  explicit  approximate  problem  satisfies  the  (first  order) 
optimality  conditions  of  the  real  problem,  that  is,  the  appro¬ 
ximate  and  real  restraint  surfaces  have  the  same  tangent  plane 
(see  Fig,  1).  As  a  result,  the  GOC  approach  should  converge 
to  a  true  (at  least  local)  minimum  weight  design.  It  will  be 
shown  in  this  section  hew  such  first  order  explicit  approxima¬ 
tions  can  be  obtained  for  various  types  of  behavior  constraints 
and  structural  models. 


4.1.  Displace  men t  Cons  traints 


The  key  idea  to  obtain  explicit  approximations  of  the  dis¬ 
placement  constraints  is  to  come  back  to  the  virtual  load  pro¬ 
cedure,  which  permits  decomposing  any  static  response  quantity 
into  the  contributions  of  each  element.  The  expression  (  b  )  can 
be  rewritten  in  the  more  general  form 
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Now,  for  a  rather  general  class  of  structural  models,  each  ele¬ 
ment  stiffness  matrix  can  be  assumed  to  have  the  following  ex¬ 
plicit  form  in  terms  of  its  design  variable  [  14  | 
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where  the  matrices  K  f  ^  are  independent  of  the 

Note  that  most  often,  at  least  one  of  the 

(2 )  .  .  1 

zero  (for  example  K.  is  zero  in  the  stiff nes 

a  flat  shell  element).  Introducing  (17)  i  n  i  n 

that  a  convenient  explicit  approximation  of  a 

traint  (5)  is  : 
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where  the  coefficients 
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are  assumed  to  be  constant  during 
dient  of  this  explicit  approximate 
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On  the  other  hand,  differentiating  (37)  and  inserting  the  result 
into  (36)  shows  that 
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Therefore  it  can  be  concluded  that  the  expressions  (38)  represent 
first  order  explicit  approximations,  in  that  they  restitute  the 
exact  values  of  the  constraints  and  their  first  partial  deriva¬ 
tives  at  the  design  point  a°  where  the  structural  analysis  is 
made  : 
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An  alternative  approach  is  to  employ  simple  linear  approxi¬ 
mations  in  the  reciprocal  design  variables  (9),  i.o.  the  first 
order  Taylor  series  expansions  (10).  Then  the  explicit  sub])  ru¬ 
ble  ni  exhibits  exactly  the  same  form  (14-16)  as  in  the  case  of 
thin  walled  structures  and  dual  methods  need  not  be  modified  as 
they  wiLl  be  in  section  5.  Of  course,  as  previously  mentioned, 
the  convergence  of  the  overall  optimization  process  might  he 
lowered,  or  even  become  unstable,  in  the  case  where  the  bending 
behavior  is  dominant  in  most  of  the  elements.  The  reader  is  re¬ 
ferred  to  sections  6.1  and  6.5  for  numerical  examples  comparing 
the  explicit  approximations  (8),  (25)  and  (38). 

4.2.  Stress  Constraints 


The  situation  is  much  more  delicate  for  stress  constraints 
than  for  disp1 aceraent  constraints.  In  contrast  to  the  case  of 
thin  walled  structures  modeled  by  bar  and  membrane  elements,  the 
stress  matrices  are  no  longer  constant  in  the  case  of  flexion- 
extension  elements.  For  illustration  consider  again  a  flat  shell 
element.  The  stresses  in  the  upper  sheet  can  be  computed  in 
terms  of  the  generalized  displacements  by 
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stress,  membrane  stress  and  bending  stress  tensors  in  member  k, 

(m) 


and  T . 


T£ “,/  ,  are  the  corresponding  stress  matrices,  which 


are  independent  of  the  design  variables.  Each  displacement  com¬ 
ponent  admits  a  first  order  explicit  approximation  of  the  form 
(see  Eq.  33) 
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Therefore,  a  natural  choice  for  the  approximation  of  a  stress 
component  could  be  as  follows  : 
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The  expression  (43)  is  still  a  first  order  approximation . 

It  is  no  longer  separable,  but  keeps  a  simple  algebraic  form 
suitable  to  specialized  algorithms  such  as  those  based  on  dual 
methods.  However,  when  considering  the  true  stress  limitation 
(e.g.  upper  limit  on  the  Von  Mises  stress),  the  explicit  cons¬ 
traint  becomes  much  more  complicate.  For  these  reasons,  and 
before  finding  something  better,  it  has  been  decided  to  folLow 
the  alternative  strategy  proposed  at  the  end  of  section  4.1, 
that  is,  to  employ  simple  linear  approximations  of  the  form  (10). 
To  compute  the  gradient  of  the  stress  constraint  with  respect  to 
the  reciprocal  variables  (9),  we  note  that,  from  the  equivalent 
Von  Mises  stress 
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On  the  other  hand,  from 
any  stress  component 


(44)  it  follows  that  the 
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y  xy 
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requency  Constraints 
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Constraints  on  natural,  frequencies  usually  consist  in  imposing 
1  owe  r  limits 
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They  are  directly  written  in  terms  of  the  squares  of  the  Ire- 
quencies,  because  these  quantities  naturally  appear  in  the 
eigenproblem  characterizing  the  structural  modal  analysis  : 


Kq.  -  wt  M  q.  =0 
J  J  J 


(  ) 


in  this  equation  K  and  M  represent  the  stiffness  and  mass  matri¬ 
ces,  and  (q.,  j  =  l,m)  are  the  modal  displacements,  i.o.,  the 

J  ...  2 

eigenvectors  solution  of  (50),  associated  with  eigenvalues  ,j. 
The  structural  mass  matrix  has  a  linear  form  in  terms  of  the 
design  variables  : 


n  n 

M  *  .  E  .  M .  +  M  =  .  £  .  a .  M .  +  M 

1=1  1  c  i  =  l  i  1  c 


(5  1) 


where  and  are  independent  of  the  design  variables. 

denotes  the  mass  matrix  of  the  ith  element  when  a.  =  1.  M  re- 

i  c 

presents  the  contribution  of  the  non-structural  masses,  such  as 
equipments,  fuel,  etc...  It  is  well  known  that  the  first  deri¬ 
vatives  of  the  frequencies  with  respect  to  the  design  variables 
are  given  by  [  see  for  example  Ref.  [  2  ],  section  4.2  )  : 
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where  m . 

J 


is  the  generalized  mass  of  the  j th  mode  : 


T 

m  .  =  q  .  M  q  . 

J  J  J 


(5  1) 


The  way  to  derive  first  order  explicit  approximations  ol  the 
frequency  constraints  is  less  apparent  than  for  displacement  cons¬ 
traints.  In  this  report,  guided  by  the  work  done  in  Ref  1  2  | 
(section  4.2),  the  following  de compos i t i on  of  the  eigenvalues  in 
terms  of  the  stiffness  and  mass  contributions  of  each  element  will 
be  used  : 
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represents  the  contribution  ul  the  non  structural  masses  to  the 


generalized  mass  in  .  (see  h  q  . 
p 1 i c i t  definitions  (37)  oi  t 

of  the  mass  matrices  M  .  ,  the 
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o  f  the  frequency  constraints 
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take  the  form  (  >  8 )  ,  wit  h 
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where  6.  -  i  only  if 
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The  coefficients  c  f  . 
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!i  d  is  0  o  the  r  w  i  s  e  . 

and  the 

mo  dified  limits  u 

P  = 


are 


1J  J 
their  v a  Lues  at  the  current  design  point.  lust  as  fo 

placement  constraints,  it  is  easily  verified  that  the 

(  J  8  )  remain  first  order  explicit  approximations  satis 


4.3) 


As  a  matter  of  fact,  they  can  be  interpreted  as 


Taylor  series  expansions  in  terms  of  ]  /  x  .  ,  )/x7,  1  /  x . 

as  independent  variables. 


(  3  h  ) 


1  ,  d  ,  3 


(  3  7  ) 


frozen  to 
r  the  dis¬ 
ex  p  r  t*  s  s  i  on  s 
lying  (42, 
first  order 
considered 


4.4.  Buckling  Constrain t s 

Just  as  the  natural  frequencies,  the  critical  load  factors 
Xj  are  defined  through  an  e i ge np rob  loro 


K  qj  - 


q  ; 


-  0 


(.  3  8  ) 


where  S  represents  the  geometric  stiffness  matrix  and  ( q  j  , 
j  =  I  ,  m )  denote  t  ii  e  eigenvectors  solution  o  )  p  r  o  b  1  e  m  (58),  asso¬ 
ciated  with  eigenvalues  the  physical  meaning  of  the  q.'s 

is  that  of  displace  me  nts  in  the  j  t  h  buckling  node  ,  for  a  criti¬ 
cal  load  factor  X  .  .  The  buckling  constraints  consist  in  imposin g 
lower  limits  on  the  buckling  loads  : 


v  .  ?*  X  .  j  =  1  ,  m  (  3  h  ) 

.1  “I 

In  this  work,  the  following  form  o  J  t  b .  ..on:  ’  rai  nLs  will  be 

adopted 
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.  ,  ,  1  1  . 

l>:(a)  —  -  T — 7— r  0  j  =  ]  ,  in  (1,0) 

.1  Aj  Va) 

because  it  has  been  found  that  better  explicit  approximations 

are  generated  when  expanding  the  reciprocal  of  the  buckling 

loads  rather  then  the  A  .  !  s  themselves. 

J 

The  stiffness  matrix  K  has  the  form  (23).  The  geometric 
stiffness  matrix  is  related  to  the  initial  stress  state  in  the 
elements  and  therefore  it  depends  implicitly  on  all  the  design 
variables  : 

n 

S  =  is j  Si  (a)  (6  1) 

It  is  worth  recalling  that  the  matrices  are  independent  o f 

the  design  variables  for  a  statically  determinate  structure. 


As  explained  in  section  3.2  of  Ref  [1 
ti vos  of  the  buckling  loads  are  given  by 


the  first  do  r i va- 


J  \  . 

_ l  _ 

3  a  . 


3  S 
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(62) 


111  opposition  with  the  static  and  dynamic  cases  previously  dis¬ 
cussed,  the  derivatives  appearing  in  (62)  are  not  directly  a  - 
vailable,  because  the  elements  of  the  geometrical  stiffness  ma¬ 
trix  are  function  of  the  stresses  acting  in  the  prebukling  state 

However,  by  assuming  that  the  terms  are  negligible,  the  gra- 

a  i 

d i e  n l s  (62)  become  easily  computable.  This  assumption,  which  is 
typical  of  optimality  criteria  approaches  for  static  constraints 
amounts  to  not  taking  into  account  the  effects  of  structural 
redundancy  : 
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(63) 


In  this  re port, the  following  decomposition  of  the  reciprocal 
buckling  loads  in  terms  of  the  contribution  of  each  element  will 
be  used  : 
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Subs  ti  tut  iiq’  t  h  e  explicit  d  e  1  i  n  i  L  i  i  -  n  (  3  ?  the  stillness  m  a  - 

trices  k  .  into  (  b  A  )  yields  the  high  q  ua  1  i  t  y  vxpli.it  a  \>  p  r  u  i  \u  a  - 
tions  of  the  buckling  load  constraints  mu  -  )  under  t  h  v  im  n  (  3  b  ) 
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Again,  the  coefficients  c  f  ?  ^ 


r  i  lied  that  t  he  e  x  p  r  n  s  s  i  o  n  s  (  \  8  )  a  re  t  i  r  s  t 
x  i  m  a  t  i  o  n  s  of  the  constraints  (  n  (J ;  ,  s  a  t  i  . 
that  the  assumption  (hi)  can  he  i n  t redo. 
( 6 1 )  (statical  dele  r  m i n  a  c  y ) . 

The  reason  why  the  buck  liny,  c..nst.  ta 
panding  the  reciprocal  critical  loads 
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e  a  t 

l  n  i 

t  s  a  .  art 
j 

n  po 

i  n  t 

.  lust  as 

n  l  s  , 

i  l 

i  s  e  a  si  1  v  v  e  - 

v  !  e 

xi»l 

i  c  i  t  a  p  p  t  <  >  - 

n  t 

*  *  » 

*> )  s  pro  v  i  de  d 

i  u  t  o 

Lh 

e  g  r  a  d  i  o  n  t 

a  re 

t  r 

e  a  Led  by  ox- 

i  ,  ■ 

•  q  • 

<)/*  ) 

i  s  that,  t  o  r 

v  3F  ) 

t  a 

k  e  s  the  1  o  r  in 

{  a  a  t 

a  I 

o  n  j’  t  h  e  s  e  a  1  i  n 

d  i  a 

s  e 

ction  J  o i 

the 

case  <>f  s  t  r  e  s  s 

v  i  t 

me 

a  n  s  that  t  h e 

d  b  v 

a 
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face  h  j  (  a  )  =  0  at  its  point  of  i  u  t  e  r  s  e  e  t  i  o  n  with  the  scaling, 

line  (see  Fig.  3  of  Kef.  12  1  )  . 


5.  SOLUT  LON  OF  Til  K  EXPLICIT  PROBLEM 


From  the  foregoing  developments,  it  appears  that,  at  each 
stage  of  the  optimization  process,  the  following  ;;  ub  p  r  o|»  1  e  n*  mu:.! 
be  solved 


minimize 


subject  to 
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i£i 
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(68) 


(Ml) 


(70) 


Unless  the  summation  on  p  in  (69)  is  restricted  to  a  single 
value  of  p  (i.e.  approximation  of  the  form  (25)  with  p  =  1,  2 
or  3),  it  is  no  longer  possible  to  find  intermediate  variables 
|  i .c .  x^  given  in  (24)  ]  in  terms  of  which  the  explicit  cons¬ 

traints  (69)  would  be  Linear.  Therefore,  the  primal  solution 
of  problem  (68-70)  is  more  difficult  to  achieve  if  a  gradient 
projection  type  of  algorithm  is  employed  as  in  the  mixed  method 
developed  in  a  previous  work  for  thin  walled  structures  |  l,  l  !>  |  . 
The  expressions  (69)  are  still  explicit  and  they  continue  to  ex¬ 
hibit  a  simple  algebraic  form.  Consequently  a  general  purpose  op- 
t.  i  m  i  *  a  t  i  o  n  al^ori  Llini  such  an  NKWSUMT  |  16  |  c  o  u  l  d  e  a  n  i  I  v  b  c  a  d  a  p  t  e  <  1  to 

take  the  constraints  (69)  into  account.  However,  because  they 
are  stiLl  additive ly  separable,  resorting  to  dual  methods  remains 
probably  the  best  strategy,  just  as  in  the  case  of  thin  walled 
structures  [1,  2,  6,  8  ] . 


The  minimization  problem  (68-70)  can  be  solved  efficiently 
as  an  auxiliary  maximization  problem  in  the  m  1  a  g  r  a  n  g i a n  multi¬ 
pliers  r'j  associated  with  the  explicit  behavior  constraints  ((>9). 
This  dual  problem  reads  as  follow  [8  ]  : 
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where  g  j  (  r  )  denote 
which  are  equal  to 


the  components  of  the  dual  function  gradient, 
the  values  of  the  primal  constraints  : 
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Setting  to  zero  the 
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braic  equations 
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h  e  i  a  e  constraint  (  7  ’> ) 

!  s ovoral  values  s  a  t i s  f v 


Note  however  that  when  s  W  v i n g  '  r 
must,  be  taken  into  accounl  and  'hat 
equation  (78),  the  one  that  u  l  n  i  ni  i  /. t  t  the  t  u  action  (7b)  must  b  e 
retained.  Standard  techniques  can  b  i ■  employed  to  deal  with  sue! 
a  simple  one- d imensionai  problem  ( s e e  c . g .  |  i 7  1 ) . 


[n  many  cases,  as  previously  meniieiud,  at  least  one  ot  the 
terms  in  the  summation  ■  n  p  in  (  7  e  \  a  a. I  (78)  disappears,  and 
it  is  possible  to  treat  the  problem  a  n i  y  t  i  c  a  1  1.  y  .  As  an  illus¬ 
tration,  consider  again  the  case  o  1  a  I  1  a  t  s  h  e  1  1  element*  when* 
the  term  i n  p  -  2  is  missing.  Problem  (7b)  exhibits  then  the 
form  (omitting  the  index  i)  : 
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or,  sotting 


(81  ) 

(82) 


The  constant  £  is  known  to  b » *  p  * '  s  j.  t  i  e  ,  a  n  d  ovt  ry  t  h  i  ng  d  e  ]>  e  n  d  s 
thus  on  the  sign  of  the  constants  ^  ‘  ^  and  C  ^  ^  ,  which  repre* 
sent  the  contribution  o f  the  membrane  and  the  bending,  respec¬ 
tively.  In  general,  equation  (82)  admits  two  solutions,  and  the 
one  that  renders  mini  m  :»  i  Lb  e  t  u  n  e  t  i  o  n  (  /  °  )  must  be  adopted. 


I’  r  o  in  the  t  o  r  e  g  o  i  a  g  d  e  v  e  1  o  o  me  n  t  s  ,  it  a  p  p  e  a  r  s  that  the  dual 
function  (71)  can  be  c  o  n  r.  i  d  e  l  e  d  as  i  function  of  the  dual  vari¬ 
ables  only.  T  h  e  d  u  a  1  p  r  c-  b  i  e  m  slat  *;  d  i  n  (71,  72)  exhibits  a  n 

attractive  feature,  namely,  it  is  a  q u a s i - un c o n s t r a i ne d  problem, 
because  taking  care  of  tin-  u v > . \ u c  g a t  i v i  r y  constraints  (72)  on  the 
dual  variables  is  s  t  i  a  i  a;  t  i  rv  *  <!  .  hi  \  n  rone  iiomuabi  t  i  vc  dual 
v  a  r  i.  a  b  1  o  s  ,  t  h  e  c  o  r  i  e  s  p  o  n  a  1  nr  ; » i  i  .>  a  i  * .  .  ■  \  \  * .  J  e .>  a  r  e  c  o  nip  u  t  e  d  by  sol¬ 
ving  {  7  b  )  o  v  t  7  8  >  a  ud  t  h  e  :>  r  i  m..  J  «  n  .  i  i  a  i  u  L  s  a  r  c  e  v  a  .1  u  a  t  e  d  l>  y  u  - 
sing  (  7  J  )  .  Tile  dual  function  (.:i  .md  its  gradient  (73)  arc  then 
directly  known  and  a  feasible  a. .cent  direction  can  therefore  be 
de  to  r  rn  i  n  e  d  .  In  a  second  ordei  .'ilpcri  ihm,  the  lu'ssi  an  matrix  of  the 


dual  f  unc  tion  has  t  o  !>  e  t,-  o  mpu  Led,: 
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and  by  differentiating  the  definition  (78)  of  a.(r)  with  respect 

to  r.  ,  it  follows  that  . 
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On  the  other  hand,  it  is  easi  lv  see:.  from  (77)  that 
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where  the  summation  on  i  is  restricted  to  the  n  active  primal 
variables,  that  is,  thoy t  that  are  not  I  i x e d  to  a  lower  or  an 


upper  limit  (of  course 


=  0  i  o  r  a  p  a  s  s  i  v  e  variable)*  K  n  o  w  i  n  g 


the  gradient  (73)  and  the  hessian  matrix  (87)  furnishes  the 
Newton  search  direction 
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where  i  is  the  stop  1  <*  r:  g  h  i  a 
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prevent  one  of  the  d  u  a  1  vari  .ib  !  e  t  r  o  m  b  e  coming  negative  |  b  ] 


The  second  order  d  u  «i  1  uptinji  zi  r  implemented  in  S  A  M  C  K  F  |  1  H  | 
has  been  e  s  o  e  c i a  1  1  y  devised  s  c »  L  h  a i  if  ;  « ■  e  k  s  the  maximum  o  t  t  he 
dual  function  by  operating  i  n  a  s  o «;  uo  iu  »■  of  dual  subspaces  with 
gradually  i  n  c  re  as  i  n  g  d  i  me  v  s  i  <.  n  .  In  t  h  i  s  way,  the  effective  di¬ 
mensionality  o  i  the  in  a  i  m  i  /  a  t  i  < »  n  j »  r  o  b  1  r  t.  never  exceeds  the  num¬ 
ber  of  active  behavj  rr  o  n  s  t  i  u  i  n  t  s  ,  w  lii  c  h  correspond  to  non¬ 
zero  dual  variables.  Past  experience  with  thin  walled  structures 
indicate  that  this  number  i  . ;  re  J  at i vc 1 y  small  in  practice,  which 
explains  the  remarkable  e  1  I  i  i  i  e  n  c  v  of  the  dual  method  |  6  J  or 
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Subsequently  the  algorithm  us  i:u'J  1  {  i  e  J  ov  int  reducing  a  simpli¬ 
fied  line  search  pro  c  e.  dure  wtu-  u  uni  L  :•  t  »•  l  e  n  g  t  h  is  a  d  u  p  t  e  d 
most  of  the  time.  A  det  a  i  led  •!  t-  a  c  r  .  p  t  i  o  n  o  i  t  h  i  s  algorithm 
can  be  found  in  K  e  t  |  t <  |  .  Hu  i . »  :  c  s  l  version,  including  the 
treatment  of  explicit  constraint."  ■>  t  the  i  on;  (69)  is  very  si¬ 
milar.  The  only  difference^  a  r .  i  !w  w*»'*  of  computing  the  primal 

variables  in  terms  ot  tut  dual  /a;  /  jbhs  (see  Kqa  (70,78))  and 

tilt?  formula  employed  f  v»r  t  i  i «.  ii  -  .  .n,  r.i t  i  *  x  I  see  L  q  (  6  /  )  J  .  It 

should  however  b  e  r  e  c  e  g  n  L  z  e  d  that  . .  •  \ <  ■  r  i  c  a  1  difficulties  might 
occur  due  to  the  lack  i  >  1  c  o  v  e  .%  ■  t  y  t  t  h  «.•  e  x  p  1  i.  c  i  t  constraints 
(69).  indeed  it  is  possible  that  the  solution  of  hq  (76)  lor  a 
primal  variable  suddenly  jumps  )  rocs  a  fixed  lower  or  upper  bound 
to  a  free  value  when  the  dual  variables  are  slightly  modified. 
This  phenomenon,  which  means  that  the  dual,  function  first  deri¬ 
vatives  are  discontinuous  along,  some  surfaces,  needs  more  inves¬ 
tigation. 


6 ,  NUMERICAL  APPLICATIONS 


All  the  examples  presented  involve  rather  sophisticated 
flat  shell  elements  that  are  characterized  by  a  displacement 
field  cubic  in  extension  and  quintic  in  flexion  (hybrid  qua¬ 
drangular  flat  shell)  I  18  ] . 

6.1.  Cantilever  Beam  with  End  Moment 


The  first  example  is  concerned  with  a  cantilever  beam  loaded 
with  a  concentrated  moment  at  its  free  end  |  see  Fig.  2  ] .  This 
problem  has  been  previously  solved  numerically  by  PRASAD  and 
H A F T K A  |  19  ]  using  an  extended  interior  penalty  function  for¬ 
mulation  for  a  beam  having  the  following  properties  :  length  = 

1U  in,  width  -  1  in,  applied  moment  =  540  in. lb  ,  Young* s  modu¬ 
lus  ®  IQ7  psi,  Poisson’s  ratio  =  0.3  and  mass  density  =  0.3  lb/in*. 
There  exists  also  an  analytical  solution,  which  was  obtained  by 
HAUG  |  20  J  .  The  analytical  optimum  design  for  a  displacement 
limit  of  0.5  in  at  the  free  end  and  specified  allowable  stresses 
of  30,000  psi  is  as  follows  (see  Fig.  2)  : 

a(x)  =  0.30  for  0  <  x  <  2.3  (90) 

a ( x )  =  0.244  x1/4  for  2.3  <  x  <  10  (91) 

It  is  worth  noticing  that  the  beam  being  statically  determinate, 
the  redesign  relation  (31)  is  exact  (with  p  =  3)  and  reduces  to 
(90,  91).  Therefore  the  optimal  solution  should  be  generated 

after  one  structural  analysis  only. 

The  structure  is  discretized  using  20  quadrilateral  plate 
bending  finite  elements  as  indicated  in  Fig.  2.  There  is  no 
design  variable  linking  and  the  problem  involves  thus  20  inde¬ 
pendent  design  variables.  The  problem  was  first  solved  by  using 
the  cubic  expansions  (25)  with  p  *  3  and  the  dual  algorithm  des¬ 
cribed  in  section  5.  As  previously  stated,  because  the  structure 
is  statically  determinate,  the  explicit  problem  (68-70)  is  exact 
and  the  optimum  design  is  obtained  in  one  single  analysis  what¬ 
ever  may  be  the  initial  thickness  a°.  Then,  the  same  example 
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flection  is  limited  to  U . 0  J  m  .  ,  \ 
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p  =  2  in  the  redesign  relations  U 1 ) ,  almost  50  iterations  were 
r  eq  u  i  red. 
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li,  J,  Simply  Supported  Square  Plate  with  Frequency  Constraint 
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shown  in  Fig.  4 . 


now  focused  on  the  simply  supported 
This  example  is  taken  from  Ref  (21 
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|  and  it.  is 
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30  . 

concerned  with  the  minimum  weight  design  of  the  plate  subject 
to  a  natural  frequency  constraint.  The  dimensions  of  the  plate 
are  10  in  by  10  in  and  its  material  properties  are  Young* s  mo¬ 
dulus  =  30  x  1  0  ^  lb/in^  and  Poisson* s  ratio  =  0.3.  In  Ref  |  2  1  I 
i't  is  stated  that  the  uniform  plate  with  pw  =  1400  is  taken  as 
the  initial  estimate  to  the  optimization  problem,  where  p  is  the 
mass  density  and  uj,the  minimum  fundamental  frequency.  Therefore, 

in  the  present  study,  assuming  steel  material  with  miss  density 
3 

p  =  0.283  lb/in  ,  the  minimum  frequency  was  chosen  as  /  2  n  - 

11.2  Hz.  Iteration  history  data  are  presented  in  Table  4  for 
three  cases  differing  by  the  minimum  thickness  constraint  (a  =  0.1 
0.05  and  0.001  in).  The  initial  design  in  cases  l  and  3  corres¬ 
ponds  to  a  uniform  thickness  a°  «  0.2  in,  while  in  case  2,  a° 
equals  0.12  in.  Final  designs  are  illustrated  in  Fig.  4.  By 
symmetry  only  a  quadrant  of  the  plate  has  to  be  analyzed  and 
designed.  It  can  be  observed  that  the  design  obtained  in  Kef 
121  ]  is  different  from  the  design  achieved  in  this  study  (case  I) 

However  both  designs  have  about  the  same  weight  (0.765  lb  and 
0.752  lb).  It  is  also  worth  pointing  out  that  the  design  of  pla¬ 
tes  in  bending  is  an  analytically  difficult  problem,  leading  to 
multiple  local  optima  and  suggesting  that  the  optimal  plate  should 
be  made  up  of  an  infinite  number  of  stiffeners  j  22  ]  . 

6.4.  I-beam  Structure 

The  third  example  involves  the  I-beam  structure  schematized 
in  Fig.  5.  The  problem  consists  in  minimizing  the  weight  ol  the 
beam  while  imposing  lower  bounds  on  the  frequencies  of  the  three 
first  eigenmodes  :  flange  flexion,  torsion  and  web  flexion.  De¬ 
tailed  data  can  be  found  in  Ref.  [2  ].  In  a  first  optimization 
exercise,  a  pure  membrane  model  was  employed.  It  involves  35 
second  degree  displacement  elements,  including  10  fictitious 
diaphragms  (without  masses).  These  dummy  members  are  introduced 
to  obtain  a  satisfactory  representation  of  the  torsional  mode. 

Only  5  analyses  are  sufficient  to  generate  an  optimum  design  for 
this  membrane  model. 


However,  when  this  final  design  was  analyzed  by  using  a  more 
accurate  model  made  up  of  flat  shell  elements,  L  h  o  torsional  fre¬ 
quency  ( in  o  d  e  2)  was  seen  to  be  violated  h y  10  7.  .  Therefore  tin 
problem  was  again  solved  with  this  now  model,  by  ‘  <  sort  im»  f  <  *  the 
theory  proposed  in  section  4.3  and  to  the  dual  optimizer  descri¬ 
bed  in  section  5  .  Iteration  history  data  are  illustrated  in 
Fig.  5  and  the  final  designs  are  given  in  Table  5  for  both  finite 
element  models  of  the  I-beam.  It  can  be  seen  that  the  use  of 
flat  shell  elements,  although  yielding  slower  convergence,  gives 
satisfactory  results. 

6.5.  U-beam  Structure 


In  an  attempt  to  consider  a  case  where  both  flexion  and  ex¬ 
tension  loadings  play  an  equally  important  role,  the  U-beam 
structure  depicted  in  Fig.  6  was  optimized.  The  2000  kg  load 
acting  at  the  tip  produces  torsion  of  the  beam.  As  a  result  the 
upper  flange  behaves  mainly  in  flexion,  the  web,  in  extension, 
and  the  lower  flange,  both  in  flexion  and  extension.  For  sim¬ 
plicity,  only  displacement  constraints  are  considered  :  the  tip 
deflection  is  limited  to  0.1  m,  while  the  relative  lateral  dis- 
placement  is  limited  to  0.005  m  (see  Fig.  6).  Three  different 
methods  are  employed,  which  ditfet  by  the  explicit  approximations 
used  for  tlie  displacement  constraints  : 


vase  1  :  full  e  x  p  a  n  s  i  o  n  (  3  ;  )  ,  yielding  problem  (  b  8  -  7  0  )  with 

p  =  l  and  3  o  n  1  y  (.the  te  r  rr  in  p  =  2  being  zero) 

case  2  :  linear  expansion  { l  >  with  p  -  1 ) ,  yielding  problem 

(14-16) 

case  3  :  cubic  expansion  (25  with  p  =  3 ) . 

1  h c  iteration  history  data  arc  given  in  Table  6  and  the  final 
designs  in  T  able  7. 


L  t  is  interesting  to  notice  that  t h e  three  methods  give 
similarly  good  results,  which  suggests  that  the  use  o  f  simple 


3  . 


linear  Taylor  series  expansion  in  terms  of  the  reciprocal  design 
variables  might  be  a  good  strategy  for  most  of  the  structural 
optimization  problems.  This  idea  was  previously  stated  by  AUSTIN 
[23  ],  but  not  believed  valid  by  many  people  when  he  published 
his  paper  ! 


ki 


7  •  C:ONCL l?S  10 N  S  AM)  K hi  COMMAN  1)  AT  1  UN  S 

It  is  now  widely  recognized  that  a  powerful  and  rather  ge¬ 
neral  approach  to  structural  o;u  imi  zai.  i*'n  is  achieved  bv  re  pla¬ 
cing  the  original  problem  with  a  sequence  of  explicit  approxi¬ 
mate  problems.  This  approach  was  initially  conceived  for  thin 
walled  structures  modelled  by  bar  and  membrane  finite  elements, 
as  well  in  the  context  of  optimality  criteria  techniques  as  in 
the  framework  of  mathematical  programming  methods  using  approxi¬ 
mation  concepts.  It  has  been  extended  in  this  work  to  deal  with 
structural  systems  made  up  of  beam,  plat  e  and  flat  shell  elements, 
with  behavior  constraints  placed  on  displacements,  stresses,  na¬ 
tural  frequencies  and  critical  buckling  loads.  The  method  pre¬ 
sented  uses  a  second  order  dual  algorithm  to  solve  each  explicit 
subproblem.  The  convergence  properties  are  independent  of  the 
number  of  design  variables,  which  is  typical  of  optimality  cri¬ 
teria  types  of  approach  as  well  as  linearization  techniques  in 
mathematical  programming.  As  a  result  large  structural  systems 
can  be  treated  at  the  expense  of  a  few  finite  element  analyses. 

Several  difficult  points  remain  to  be  clarified.  First  it 
is  not  well  known  whether  the  explicit  problem  can  be  solved 
conveniently  in  any  case,  because  its  lack  of  convexity  might 
lead  to  discontinuity  in  the  dual  function  gradient.  Secondly, 
the  use  of  simple  linear  expansions  with  respect  to  the  recipro¬ 
cals  of  the  element  transverse  sizes,  gives  rise  to  satis  factory 
results  in  many  cases.  This  suggests  that  more  complicate  expli¬ 
cit  approximations  such  as  those  proposed  in  this  work  mi  g  lit  not 
be  necessary.  Finally,  there  remains  the  question  of  the  stress 
constraints,  for  which  it  is  difficult  to  generate  adequate  first 
order  explicit  approximations.  Also  stress-ratioing  algorithms 
should  be  considered  for  extension- flexion  elements,  because  they 
correspond  to  much  more  inexpensive  zero  order  approximations. 


REFERENCES 


34  . 


|  1  ]  C.  FLEURY  and  G.  SANDER,  "Structural  Optimisation  by  Finite 

Elements",  AFOSK  Final  Scientific  Report,  Grant  77-  MIR 
LTAS  Report  SA-58,  University  of  Liege,  1978 

[  2  ]  C.  FLEURY  and  G.  SANDER,  "Generalized  Optimality  Criteria 

for  Frequency  Constraints,  Buckling  Constraints  and  Bending 
Elements" ,  AFOSR-TR-80-0 107,  1979 

[3  ]  L.  BERKE  and  N.S.  KHOT,  "Use  of  Optimality  Criteria  Methods 
for  Large  Scale  Systems",  AGARD-LS-70,  1974,  pp.  1-29 

(4  1  L.A.  SCHMIT  and  H.  MI  UR A  ,  "A  New  Structural  Analysis/ 

Synthesis  Capability  -  ACCESS  1",  A 1 AA  J . ,  Vol.  14,  No. 5, 
1976,  pp .  661-671 

IS  1  L.A.  SCHMIT  and  H.  MIURA,  "An  Advanced  Structural  Analysis/ 
Synthesis  Capability  -  ACCESS  2",  Int ,  J.  Num.  Me th .  Engng. 
Vol.  12,  No. 2,  1973,  pp,  353-377 

l  6  ]  C.  FLEURY  and  L.A.  SCHMIT,  "Dual  Methods  and  Approximation 
Concepts  in  Structural  Synthesis",  NASA  CR-J226,  1980 

[7  ]  R.L.  FOX,  "Constraint  Surface  Normals  for  Structural  Synthe 
sis  Techniques",  A I A  A  J ,  ,  Vol,  3,  No. 8,  1965,  pp.  1517-1518 

(  8  ]  C.  FLEURY,  "Structural  Weight  Optimization  by  Dual  Methods 
of  Convex  Programming",  Int.  J ,  Num.  Methods  Engng.,  Vol. 14 
No.  1  2  ,  1979  ,  pp  .  1  761-  1  7  8  3 

[9  )  J.  KIUSALAAS,  "Minimum  Weight  Design  of  Structures  via  Op¬ 
timality  Criteria",  NASA  TN-D-7115,  1972 

[  10]  A . J .  MORRIS,  P.  BARTHOLOMEW  and  J.  DENNIS,  "A  Computer  Base 
System  for  Structural  Design,  Analysis  and  Optimization", 
AGARD  CP-280,  1980,  paper  20 


[  II  ]  C  .  PE’ilAC  and  G  .  Idit  INA,  "  1.  i  C  me  u  L  s  Finis  it  Optimisation 
des  S  triK'tures"  ,  A  G  A  K  i)  C  V  -  7  8  0  ,  1  a  8  0  ,  p  a  p  e  r  2  3 

1  12  1  J.L*  A  KM  AND  and  B.  1.0!'  IKR  ,  "  o  p  t  i  ma  1  I)e  :•  i  v.  n  nl  Ben  din/ 

Elements",  1  at  .  2L.  Num.  M  u  th  ods  hngrp].  ,  V  o  l  .  I  3  ,  N  o  .  2  , 

1978,  pp.  3/3-384 

|  1  1  1  W  .  P  R  A  G  E  K  ,  "Optimal  Design  ot  Static  a  1  1  v  l)o  to  r  m  L  n  a  t  o  B  o  a  in  s 

1  or  Given  Def  It-  tiun"  ,  In  t_. _ J  .  M  e  c  h  .  S  c  i  .  ,  ini  ,  1  >  ,  j  a  7  ]  f 

p.  89  3 

|  14  ]  J  .  K  1  l’  S  AL  A  AS  and  K  .  C  .  J  .  5  H  AW,  "An  A  1  g  o  r  i  t  h  m  lor  0  p  t  i  m  a  1 

Structural  Design  with  Frequency  Constraints"  ,  1  nt  .  A  . 

N  uni,  Me  thods  Kngng  .  ,  V  o  i  .  I  -I  ,  t\  o  .  2  ,  1  0  7  8  ,  p  p  .  2  8  3  -  2  9  > 

1  15  ]  G.  SANDER  and  C.  FLEURY,  "A  Mixed  Method  in  Structiw.il 
Optimization",  1  n  t^.  *  *  •  ^  u  111  -  M  e  t  li  .  tngng.  ,  Pol.  1  i  , 

1978,  pp.  3  85-40-* 

|  16  1  L  .  A  .  S  C  HM 1  1  and  li  .  M  I  F  R  A  ,  "Approximation  concepts  to  r 
Efficient  Structural  Svntlu  s  is",  CASA  c  K  -  2  ' ,  1  1  *  7  ♦  - 

(17  1  M.A.  WOLFE,  Numerical  Methods  lor  inconst  r.i  i  nt<  a  U  p  t  i  m  i  z  a  - 

tion  -  A  n  Introduction,  V  a  n  N  o  s  t  r  a  n  u  K  o  i  nh  o  l  d  %  W  o  k  i  n  g  h  a  m  , 

1  9  78 

1  18  |  S  A M C E F  ,  S  y  s  l  e  m  e  tl  ’  A n a  lyse  dts  Milieux  l  or,  t  i  nus  par  I,  I  emu 
F  i  n  i  s  ,  LIAS,  University  of  1,1c  >•  t  ,  B  e  1  g  i  uni 

|  19  )  B.  PRASAD  and  R.T.  HAFTKA,  "Optimal  Structural  Design  with 
Plate  Finite  Elements",  .  K  S u c  t  u  r  a  1  Pi  v  i  s  i  on  *  A S  C F.  , 

Vo  1  .  105,  No. STM,  19  79  ,  pp.  2367-2982 

(20  1  E.J.  HAUG  and  P.G.  K1RMSEK,  "Minimum  weight  Design  ol 

Beams  with  Inequality  Constraint  s  o  u  8  t  r  e  s  :>  and  I)  e  J  J  *  c  1  i  «*  n 

J.  Applied  Me  ch  a  n  i  c  s  ,  19  6  7,  p  p  .  99  9  -  1  ( 1 0  4 


E. J.  H AUG ,  K  .  C  .  PAN  and  T.D.  STREETER,  "A  Computational 
Method  for  Optimal  Structural  Design.  1.  Piecewise  Uni¬ 
form  Structures"  ,  Int.  J.  Num.  Methods  Engng.,  Vo  I.  i. 

No  .  2  ,  1972,  pp.  171-184 

N.  OLHOFF  and  KENG-TUNG-CHENG  ,  "Optimal  Design  of  Solid 
Elastic  Plates",  in  Optimization  of  Distributed  Par  a  m  c  U?  r 
Structures ,  (E.J.  HAUG  and  J.  CEA,  eds),  Sijthoff  and 

Noordhoff,  1981 

F .  AUSTIN,  "A  Rapid  Optimization  Procedure  lor  Structures 
Subjected  to  Multiple  Constraints",  Proc.  A  L  A  A  /  A  S  M  K  /  SA1' 
18th  Structures,  Structural  Dynamics  and  Materials  Confe¬ 
rence,  San  Diego,  California,  1977,  pp.  71-79 


Table  2 


Final  Designs  for  Cantilever  Beam 


E  1  erne  n  t 
No. 


Weight 

(lb) 


Thickness 

(in) 

Analytical 
Ref  [  20] 

Ref  [  19  ] 

Linear 
Expans  ion 

0 . 3 

0 . 3004 

0  .  3 

0  .  3 

0.3006 

0  .  3 

0. 3 

0.3018 

0  .  3 

0 . 3 

0 . 3056 

0 . 3 

0 . 3 

0.3109 

0.  3 

0.3142 

0.3167 

0.3167 

0.3276 

0 .323  1 

0.3287 

0. 3396 

0 . 3298 

0 . 3327 

0 . 3503 

0 . 3367 

0 . 3350 

0 .3602 

0.3437 

0.3379 

0  .  3693 

0. 3508 

0 . 3450 

0.3778 

0.3579 

0.3574 

0 .3858 

0 .3652 

0.3711 

0.3933 

0 .3726 

0.385  1 

0 . 4004 

0.3803 

0. 3989 

0.4071 

0. 3881 

0.4121 

0.4135 

0 . 3960 

0 .4246 

0.4197 

0.4040 

0.4  358 

0.4255 

0.4  119 

0.4446 

0.4312 

0.4  197 

0 .4533 

1  .0750 

l  .0453 

1  .  0768 

0  , 

.31  36 

0  , 

.  32  70 

0  , 

.  3  389 

0  , 

.34  9  7 

0  . 

.  3596 

0  . 

,  3687 

0  , 

.37  71 

0  , 

.  385  1 

0  , 

.  39  26 

0  . 

,  3996 

0  , 

.406  3 

0  . 

,4126 

0  . 

.4182 

0  ■ 

,4  225 

0  , 

,4268 

1.0797 


No  -  of 
Analyses 
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8 
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Table  4  :  Iteration  History  Data  for  Square  Plate 

with  Frequency  Constraint 


Analysis 

No. 


Case  1  :  a  =  0  .  1  in  Case  2  :  a  =  0 . 0  5  in  Case  3  :  a  =  0 . 0  0  I  in 


Weight 

(lb) 

Frequency 

(Hz) 

'  Weight 
(lb) 

! _  .  .  . 

F  req  uen 
(Hz) 

1.415 

19.6 

0.849 

11.8 

1.097 

16.0 

1  0 . 786 

11.5 

0.913 

13.6 

0. 750 

11.5 

0.812 

12.2 

0. 708 

11.7 

0.  771 

11  .  5 

0.658 

11.6 

0.757 

11.3 

0.63  1 

11.4 

0 . 753 

11.2 

0.617 

I  I  .  3 

0.752 

1  1  2 

0.605 

11.3 

0.752 

11.2 

0.596 

11.3 

0.589 

11.3 

i 

0.583 

11.3 

i 

0.5  77 

11.3 

IBM  370-  1  58 
CPU  time 
(sec) 


1.415 
1  .097 
0 . 902 
0.778 
0.696 
0.647 
0.622 
0.604 
0 . 588 
0.574 
0.559 
0.542 
0.52  1 


nge  d 


Table  5  :  Filial  Designs  for  I-beam  Structure 


thickness  ( mm) 


membrane  model 
i 

flat  s  lie  1  1  model 


Upper  Flange 


15.80 

17.30 

11.67 

6.14  3 

1.813 

16.59 

15.09 

10.08 

5 . 090 

1  .  36  3 

_ 

Web 

5.101 

3 . 602 

3 .329 

3 .2  94 

1  .  9  9  7 

2  .056 

4.173 

3.965 

3.930 

2.313 

t- 

S  upported 

Lowe  r  FI ange 

15.40 

17.01 

11.5  5 

6.074 

1.79  3 

20.51 

17.98 

i  2.21 

6.480 

1.91  () 

Table  6  :  Iteration  History  Data  for  U-Beara  Structure 


NUT** 


Table  7 


Final  Designs  tor  11  -  IS  e  a  m  Structure 


V  p  p  e  r 
F i ange 


Web 


Lower 


Thickness  ( mm ) 


Case  l 
C  a  s  o  2 
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Full  Expansion 
Linear  1  x p a  ns i n 
Cub  i  c  E  a  p  an  s  i  <  >  n 
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FIG.  2  CANTILEVER  BEAM  WITH  END  MOMENT 
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SIMPLY  SUPPORTED  SQUARE  PLATE 
WITH  DEFLECTION  CONSTRAINT 
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FIG.  5  I-BEAM  STRUCTURE 
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FIG  6  U-BEAM  STRUTURE 


